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Abstract 

In this paper we derive a general expression for the emission rate of an emitter which is placed 
at the center of a real cavity drilled in a generic random medium. We apply our formalism to the 
computation of the decay rate of a single atom which seats within a small molecule embedded in 
a continuous medium. 
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In a previous paper [3| we computed the spontaneous emission rate of a fluorescent 
particle seated at the center of a spherical cavity of radius R drilled in a continuous medium 
for the case that R satisfies R <^ X, R ^ ^, with A being the wave length of emission and 
^ the correlation length of the medium constituents. We referred to the former inequality 
as small cavity limit. For the case that both inequalities hold, the whole system cavity- 
host-medium can be treated as a continuum and a homogenous bulk propagator can be 



computed in good approximation. In we used the formal expression obtained in p| for 
a real cavity scenario in the small cavity limit. In the present paper, we will relax the 
conditions about the continuity of the system. We will still consider that the host medium 
is continuous as seen from the emitter beyond r > R. That is, R^ ^. We will refer to this 
inequality as macroscopic cavity limit. However, the freedom in the choice of the ratio A/^ 
implies that spatial dispersion in the dielectric constant is allowed. Secondarily, the lack of 
constraint in the choice of the ratio X/R implies that the emitted wave does not only see 
the cavity twice in the processes of emission and back-reaction, but also from any scattering 
event within the host medium. Propagation rather depends continuously on the cavity and 
a translation-invariant bulk propagator in the medium surrounding the emitter cannot be 
strictly computed. 

Our first goal is to derive a recurrent formula which accounts for the whole series of 
radiative corrections over the bare polarizability function of a point emitter. On the other 
hand, we would like it to be comparable with the exact formula found in ll| for a virtual 

n 

cavity scenario. Following p| we will consider a point polarizable emitter. Its polarizability 
function a{k) gets renormalized as a result of iterative self-polarization processes, 

aCk) =ao[l + ^aoP3fJ{27f } + ^ao~k'{2^± + 7||)]~^ (1) 

In the above equation, the 7-factors encode the radiative corrections over the bare polariz- 
ability The term ^aok'^^{2'y^f^} has been introduced to account for the internal reso- 
nance of the emitter in vacuum. As shown in jsl , it plays the role of a regulator of the intrinsic 
ultraviolet divergence in . That is, ?fi{2-ff^} = where tto is the corresponding real 
electrostatic polarizability and ko is the resonance wave number in vacuum. Following [6], by 
parametrizing Eq.([l]) in a Lorentzian (L) form, aL^k) = a'^kl^Jikl^^ — k"^ — iV^k^ / {ckl^g)]~^ , 
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we can identify the decay rate of the emitter in the host medium as 



pa 



--a[,fc^53{27x + 7||}k=fc.. 
-ro|jfc53{27^ + 7||}|^=,., 



where kres is a real non-negative root of the equation 

1 



(2) 



aok'n^^^ + l\\}\k=k. 



(3) 



«o ~ (^oiko/kresY IS the renormahzed electrostatic polarizability and Fq = caoA^o/Svr is the 
in- vacuum emission rate. As argued in 6|, consistency with Fermi's golden rule requires 
ao = -^j^i^, A* being the transition amplitude between two atomic levels in vacuum. 
In the case that the emitter itself is a host scatterer, the emitter occupies a virtual cavity 



and the factors 7^j| are given by 

27± = 
7|| = 



d^A; 

(fk 
(2^ 



2x±(fc ), 
pa 

X\\{k) 



G^{k) -23?{7f} 



pa 



Gn{k) 
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(4) 
(5) 



In the above formulae, x{k) is the effective susceptibility and G±,|| are the transverse and 
longitudinal components of the bulk propagator respectively. 

In the present case, the emitter seats in a real empty cavity of radius R. Therefore, it 
does not correlate to the host medium as any host particle would do. Assuming R ^ ^, 
the complete series of diagrams is that in Fig{T], where we use the diagrammatic conventions 
employed in in the macroscopic cavity limit. 

The first thing to notice is that the condition R^ C, allows for simplifications. Let g{r) 
be the two-point correlation function of the emitter with the host scatterers. 



g(r) = 1 + hc{r), 



(6) 



where hc{r) is the irreducible piece. Because the cavity is macroscopical, it distorts in 
general the homogenous distribution of host scatterers in the surrounding medium. Thus, 
hc{r) must be provided. For the moment, let us think of it as a step function which accounts 
for the cavity exclusion volume, hc{r) ~ — B(r — R). Its Fourier transform reads 
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{k) = {27i)H'^^\k) + hc{k). 



(7) 
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FIG. 1: (0) Feynman's rules, (a, 6) Diagrammatic representation of the transverse and longitudinal 
7-factors, according to Eqs. (|ll|12p (a) and Eqs. (|13|14p (6) in the macroscopic cavity limit. 

Because hc{r) has support in r < i?, it is therefore expected that its Fourier transform gets 
support for k <1/R. On the other hand, in case that spatial dispersion in e{k) be relevant, 
such a dispersion must be of the order of 1/^. Therefore, any convolution of g{\k — k'\) with 
the self-energy operator and any generic function f{k') can be approximated by 

S(A;) is a tensor made of one-particle-irreducible (IPI) multiple-scattering processes which, 
for a statistically homogeneous host medium, can be decomposed as S(A;) = T,±{k){l — ^) + 
In Eq.(IH]), the difference between k and k' is negligible in comparison to 1/^ in 
the range of momenta where hc{\k — k'\) takes nearly constant value and ^{k) is not zero. 
In turn, that implies that the correlation of the emitter with any IPI (multiple) scattering 
process can be approximated by the correlation of the emitter to any of the scatterers 
involved in such a process. In particular, things get mathematically simpler if correlations 
connect either the entrance or the exit scatterer in the corresponding IPI diagrams. This 
allows to write the series of radiative correction as two apparently different expansions, (a) 
and (6), as depicted in FigsII](a, 6). The corresponding series read 



n=0 n=0 

oo oo 



n=l n=l 



In the above series, the n*'* term of each one reads 



27 



a{n) 



(0), 



(27r)3 

d^A; 



(11) 

(12) 



and 



27^") = 2 
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b{n) 



d-'k .(n-l). 
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-x\^\k)^ik)G^;\k)^ik)xl;'ik). 



(13) 
(14) 



The recurrent formulae for the partial pseudo- susceptibilities X^± \ -the reason for this nomen- 
clature will get clear later on- read 



3u 



d-'k 



x[r^(^) 



^2^^s9i\k'-k\)[il + cos^ 
X xt'\k')^^{k')G^l\k') 
+ sin^e x^f-'\k')mk')G^p{k') 
d^k' 



(15) 



(27r) 



cos^^ 



^,9{\k'-k\) 
X xt'''ik')mk')Gf\k') 



+ sm'ex'-r'\k')j:^{k')G'l>{k') 



(16) 



for n > 1 and X^j^^^^{k) = — A;^S_|_ ||(A;) ^ for n=0. In Eqs. fll5lll6p . ||(A;') can be factored 
out of the integral as S_lj|(A;) in application of the approximation in Eq.([8]). We notice that 



xli^l(fc) relates to the cavity factors firstly defined in [ij 



C^k) 



dV e'^-''hc{r)TT{G^^\r)[I - k k]} 
^'^\hc{\k'-k\)\G^l\k') 



Ciiik) 



(27r)3 

+ G^^\k') cos^e + G\^\k') sin^e], 
dV e'''-''hc{r)TT{G^^\r)[k ® k]} 

§^shci\k'-k\) 
Gf\k') cosH + (fcOsin^^' 



(17) 



(18) 
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by X±\\{k) = '~^'^[G±\{k) + C±^\\{k)] in a homogeneous host medium. Also, Eqs. fll3|14p 
resemble the formulae of the 7-factors for the virtual cavity scenario -Eqs. (l4|5|) . That is, if 
we define the total pseudo-susceptibility as x±,\\{k) = ^n=oX±\ik) we can write 

However, the above equations are not yet quite similar to those for the virtual cavity scenario. 
In particular, the propagator in the integrand of Eqs. ( fT9]l20l) is that of free-space whereas it 
is the bulk propagator for the virtual cavity scenario. The reason being that x±^\\{k) contain 
both IPI and non-lPI processes. We can decompose x±,\\{k) in IPI and non-lPI (NIPI) 
pieces, x_L^\\{k) = xY'^^{k) + X^Y^(^) according to the following decomposition in partial 
pseudo-susceptibility functions: 

xZr^\k) = x^rj'\k)^U\ik)GTik), n > 1, 

n > 1, 



x{'i^-\k) = tlik)-xZr\k), 

= xZik), x:r^fc)=o. (21) 



Using this decomposition, we can write Eqs. (ll9|20l) in the form. 



^-^^/(Slw^^'^-*' '''' 

(fk x]r{k) . 



where 



G^,\\{k) ^ G%{k) + Gr||(^)^^^TW77:/s^Jl(^) (24) 



From this equation we can identify a T-matrix for the emission process, 

T±,\\{k) = ^%jirv^um. (25) 



XI 



Applying the same reasoning as that in 



m 



we decompose the above function in prop- 



agating (P) and non-propagating (NP) components -in absence of longitudinal propagating 
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modes. That is, 

The above equations are suitable for any real cavity satisfying i? ^ ^. A generic case 
will be treated numerically in a separate paper. Here, in order to make contact with pre- 



vious results 



0,3 



we will treat a simpler realistic scenario. That is, we will compute the 
spontaneous emission of a two-level atom which seats within a small molecule embedded 
in a continuous host medium. By small molecule we mean that the molecule radius Ri is 
such that -Ri ^ A^es- By continuous medium we mean that Ri ^ ^ being ^ the correlation 
length of the host medium constituents as in [s] -see FiglHl 

Consider first the molecule is in vacuum as depicted in Figl2](a). The molecular medium 
surrounding the emitter atom can be thought of as a continuous effective medium of permi- 
tivity ei within Rq < r < Ri. Hence, the cavity correlation function reads 

gc{r) = e{r-Ri)-e{r-Ro). (28) 

Thus, gc is totally irreducible and so are x±\\{k) and x±,\\{k)- As a matter of fact, 

xi\)(A;) = -P[C'5|-C5|], (29) 

where the superscripts denote the radii of the cavities and Eqs. fll7|18p are used. Further on, 
in the small molecule limit -Ri <^ A^es, the wave cannot resolve the molecule structure at all, 
so that the field is roughly uniform in Rq < r < Ri and we can write -see Figsl2](c, d)- 

X±,\\{k) = + X^Yrn - X%{k)/x^^Y'- (30) 

That is, we can interpret this result as if the wave visited infinite times the molecule before 
propagating throughout free space and got back. In this case -and only in this case, the 
7-factors reduce to those formulae in Eqs.( |T9]l20l) with X-L,||(^) given by Eqs. ( 129)1301) and 
X^O) = l/(ei — 1). From those equations it is easy to verify that the 7-factors can be 
expanded in powers of (ei — 1). 27^(A;) can be computed analytically in closed form. We 
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FIG. 2: (a) Schematic picture of a point emitter seated at the center of a small molecule of radius 
Ri in free space. Rq is the radius of the emitter cavity in the molecular medium. (6) Diagrammatic 
representation of the actual expansion of the 7-factors according to Eqs. (j9lllll2p . (c) Diagram- 
matic representation of the 7-factors expansion in the small molecule limit, (d) Diagrammatic 
representation of Eq. (j30p in the same approximation, (e) Diagrammatic representation of the 
summation of the above series in the small molecule approximation in free space. 

give below the leading order terms up to order (ei — 1)^, 

i r (0)t 



i r 



i / 1 1 \ Ilk 



2~k\Ri Rj 30 
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It is remarkable that, in the small molecule limit, the Lorentz-Lorentz (LL) local field factors 
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cancel out with respect to the computation of [3| for i?i ^ oo and in accordance with the 
numerical computation of j4|. Also, for a given power n of (ei — 1), the leading term in (kR) 
is of the order of 2n but for n = 1 which contains a non-propagating contribution which goes 
as ~ j^. On the contrary, the series expansion of 7||(A;) contains near field contributions at 
all order n in (ei — 1). That is, we can write the series for 7||(/c) as 



k 



27r 

n=l 

^(n) ^(n), 

kRi kRo 



{kR,y {kRo^ 

(32) 



where A^^^ and b["^ are real numbers of order 1. In particular, up to order n = 2, A^^^ = 
l,Bf^ = -1,4'^ = 1/2, = -1/2,4'^ = -2/3, 5f = 2,4'^ = 0, ^ =4/3. 

Let us introduce next the molecule of radius Ri inside a continuous medium of permitivity 
e2 as shown in Figl3l This way, the scenario is similar to that addressed in jsl in which the 
emitter seats at the center of a small real cavity drilled in a continuous medium. The 
difference here being that the cavity itself is a molecule which hosts the atomic emitter. 

n 

The series of diagrams is that of FigJHl It is analogous to that in [3] but replacing the 
factors Cj||(fc) + Gf^^^{k) with ^^^[Cf^i^) + G^J^^A;)]. We will restrict ourselves to the 
computation of the propagating decay rate, 2V^{k) = ^Y^To'^{2'y^}(k), where k is the 
resonant frequency in the medium satisfying Eq.([3]) and Tq is the decay rate of the atom in 
absolute vacuum. '^{2'y^}{k) reads, 

-2P5J{|^(e,-l)Gf (^)} I ^{G'l\k)}^, 



(27r)3 

3) 



X. J 
+4k'^{^{e,-l)C^^Ck)} J '^{GTi^2-l)G^{k)}-^ 

where G±{k) = ^ j}_f^2 is the transverse component of the bulk propagator in the continuous 
medium. For the sake of simplicity, let us assume ei,e2 are real in first approximation. In 
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FIG. 3: Diagrammatic representation of the 7-factors for the case of a point emitter seated in a 
small molecule embedded in a continuous medium of permitivity €2- The propagating transverse 
component is that in Eq.(| 



such a case, we obtain 



- i(e2 - l)Pi?i)'(4 - 61) + (61 - mRof] 



X [l + -7i;-2(e2-l)(v^-l)] 



(34) 



Where 2r^ is the formula found in yl when the finite size of both the molecule and the 



emitter cavity are discarded. It does contain the LL local field factors 



2T^ = To 



'2 + 2^^ h 



(35) 



The additional terms account for finite size corrections at leading order in kR. As argued 
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above, the expansion of x±/x^^^ i^i powers of (ei — 1) yields terms ~ (ei — at 
leading order in {kR). Therefore, in the range of frequencies where the small molecule limit 
is valid, strong deviations are expected with respect to Eq. fl34|) close to resonances of the 
permitivity of the intramolecular medium. That is, for frequencies ck such that kRi <^ 1 
but|(ei-l)(^i?i)2|>l. 

In summary, we have derived generic expressions for the 7 radiative factors which dress 
up the bare polarizability of a point emitter which seats within a macroscopic real cavity 
{R ^ C,) drilled in a random medium -Eqs. fl22|23l) together with Eq. (1241) . To that aim 
we define pseudo-susceptibility functions, x_L,\\{k), whose partial components satisfy the 
recurrent relations in Eqs.( fT5|T6l) . The formulae so obtained are analogous to those found 
m in a virtual cavity scenario. Propagating and non-propagating emission rates are 

identified -Eqs. ( 126)1271) . For the case of a single atom which seats within a small molecule 
embedded in a continuous medium an expression for the pseudo-susceptibility is found in 
closed form -Eqs. (!29|30l) - in terms of the cavity factors -Eqs. (ITTllTSl) . An expression for the 
propagating emission rate is found as a function of the radii of the molecule and the emitter 
cavity within the molecule and the permitivities of the intramolecular medium and the host 
medium -Eq. . 
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